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Abstract 

In ||] a new bosonization procedure has been illustrated, which allows to 
express a fermionic gaussian system in terms of commuting variables at the 
price of introducing an extra dimension. The Fermi-Bose duality principle 
established in this way has many potential applications also outside the context 
of gauge field theories in which it has been developed. In this work we present an 
application to the problem of averaging the correlation functions with respect 
to random potentials in disordered systems and similar problems. 
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1 Introduction 



Both in statistical mechanics and in quantum mechanics there are several sit- 
uations in which one has to average the correlations functions of a physical 
system with respect to disorder fields The averaging procedure is com- 
plicated by the fact that the dependence of the correlation functions on the 
disorder is not explicitly known. If we restrict ourselves to systems which ad- 
mit a representation of the Green functions in terms of gaussian fields, one has 
two powerful tools at disposal to solve this problem, the replica method [|| and 
the supersymmetric method The replica method can be applied to more 
general cases, but has the disadvantage that its mathematical consistency has 
not yet been proved, so that it has been subjected sometimes to some critiques 
concerning its validity in nonperturbative calculations On the other side, 
the introduction of supersymmetric fields requires a fermionization of the sys- 
tem, in which the passage from bosonic to fermionic degrees of freedom is often 
obscure. For this reason, difficulties arise for instance in systems with sponta- 
neous symmetry breaking, because it is not easy to interprete the symmetry 
breaking in terms of the resulting fermionic theory ||. 

In this letter we propose an alternative to the above two methods. Our 
approach is very similar in spirit to the supersymmetric one, but it involves 
only bosonic fields. Instead of fermionization, we exploit the Fermi-Bose dual- 
ity principle introduced by A. A. Slavnov in the context of lattice gauge field 
theories in order to provide an expression of fermion determinants which does 
not contain anticommuting fields Jl|. Successively, this principle has been also 
applied in to rewrite the Faddeev-Popov determinants without anticommut- 
ing ghosts. Since in our case all operators are hermitian, it has been possible 
to simplify the original procedure, which otherwise would have lead to field 
theories with derivatives of the fourth order in the action. Another relevant 
change with respect to |7| is the choice of boundary conditions of the auxil- 
iary bosonic fields. Here boundary conditions are dictated by the compatibility 
requirement with the regularization needed in the path integral approach to 
quantum mechanics to guarantee convergence. 

The material presented in this work is organized as follows. In Section |2| 
the problem of averaging over the disorder fields is briefly discussed. In Section 
|3| our alternative method based on bosonization is presented. The conclusions 
are drawn in Section |||. 

2 The Averaging Problem In Disordered Sys- 
tems 

Let H be a local and hermitian Hamiltonian describing a system with D degrees 
of freedom q = (g 1; . . . ,qn)- Further, we suppose that H depends on a set of 
random potentials v?(q) = (y?i(q), (p2(q), . . . ) with a given distribution P(0). 
Hamiltonians of this kind are widely applied in quantum mechanics 0. For 
instance, choosing <p = (tpi), one obtains the energy operator of a disordered 
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system 

H = H + ip! (1) 

consisting of a fixed Hamiltonian Hq and a random perturbation y?i(q). If for 
example <*p\ is a source of Gaussian noise, the general form of P(<pi) is given 
by: 



p(^) = ex P |- J dV D gVi(q)^(q,q>i(q': 



) (2) 



Analogously, with the same formalism it is possible to discuss the motion of 
n— dimensional particles immersed in an electromagnetic field with components 
Ai, i = 1, . . . , n, putting if = (A\, . . . , An) and taking as "distribution" 

P(A 1} ... ,A D ) =exv{-iS QED } (3) 

where 

Sqed = ^ J d D xFl (4) 

is the usual action of quantum electrodynamics in n dimensions and Fij = 
diAj djAi. 

In the following, we denote the quantum average and the average over dis- 
order fields with the symbols ( } and ( )^ respectively. With this notation, 
the advanced and retarded Green functions of the Hamiltonian H are given by: 



G E {q_,c{;0) = lim ( q 



1 



E±ie-H 



q' ) (5) 



where E + ie is a complex parameter with an arbitrary small imaginary part ie. 
Relevant information about the system may be obtained computing averages 
over (p of products of the above Green functions: 

(G±(q,q , ;^)G±(q,q , ;^)---)^= j V0P(<p)G%(q, q'; 0)G% (q, q'; 0) ■ ■ ■ 

(6) 

To this purpose, it is often convenient to use a representation of G^(q, q'; V) 
in terms of complex scalar fields <f>, (j>: 

G±(q, q'; 0)) = J P0P#(q)0(q')exp l±i J d D q0(E ±ie- H)A (7) 

Z± represents the partition function of the field theory: 

Z± = J V(j)V4>exp |±i J d D qj)(E ±ie- H)A (8) 

It is easy to realize that, even in the simple case of a single noise source with 
gaussian distribution as in Eq. (pi), it is difficult to integrate over the random 
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potentials in the right hand side of Eq. (||) due to the presence of the factor Z^ 1 
in the definition of G ± (q, q', E\(p) (see Eq. ([?])). In fact, the partition function 
Z± is a functional depending on (p in a complicated way. In the next Section, 
it will be shown how it is possible to perform the average with respect to the 
random potentials without introducing replica fields or fermionic degrees of 
freedom. 



3 The Bosonization Method 

First of all, we note that the class of problems under investigation has a gaussian 
nature, as it is shown by the field theory representation of Eq. (0), in which only 
gaussian fields are involved. Due to this fact, for our aims it will be sufficient to 
consider only the average of a single Green function. Let us study for instance 
the following average: 

(G-(q,q'; J V0P(<p)G-(d,q!;<p) (9) 

Now it will be convenient to interpret the factor Z^ 1 in Eq. (Q) as the functional 
determinant of the operator E — H: 

Z^ 1 = det(£ ±ie-H) (10) 

To express the determinant appearing in the right hand side of Eq. ([h]), we 
apply the Fermi-Bose duality principle proposed in Q. To this purpose, we 
introduce a fictitious time r such that 

—T < t < T (11) 

and two sets of auxiliary fields c n (q, r), c n (q, r) and Xn(q), n — 1,2. The field 
c n is the hermitian conjugate of c n 

Cn = (Cn) ] (12) 

while Xn is an hermitian scalar field, i.e. (x-n)^ = Xn- The fields c n and c n 
satisfy the boundary conditions: 

c n (q,-T) = B n (d) c n (q,T) = £?t(q) (13) 

where -B ra (q) is an arbitrary function of q. 

We are now ready to prove the following formula: 

ZZ 1 = det(E + ie - H) = lim Z cl Z c2 (14) 

where 

Z c , n = I Vc n Vc n V X ne iSn (15) 
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and 



S n = J T dT J d D q 

i dc n 



2 dr 



(E-H) + 2iec n c n + x( c n + Cf> 



(16) 



The proof goes as follows. Since the operator E — H is hermitian, it supports 
a complete system of orthonormal eigenfunctions ip a with eigenvalues X a : 

(E — H)if) a (q) = X a ip a (q) (17) 

Thus, it is possible to expand the fields c n , c n and Xn in terms of the ip a : 

Cn(q,r) = £c"(rM»(q) (18) 

a 

c n (q,r) = ^c^(r)^(q) (19) 

a 

Xn(q) = 5>^«(<0 (20) 

a 

Here c"(t),(^(t) depend only on the pseudo-time r, while the Xn are constant 
coefficients. To these equations, one should add the expansions of the fields B n 
and Bn which express the boundary conditions: 



B n ( q ) = KM<\) 4(q) = E Ki>M) 

a a 

Substituting Eqs. (l||-2C) in the action (16) and remembering that: 

j d D qij a (q)ipp(q) = 5 a(3 

one obtains for Z cn : 



(21) 



(22) 



(23) 



with 



(It 



T 



2 dr 



+ A a c! 



i dc° 



+ ( ~~ ^ aC " ' C ™ + ^ eC n C n + Xn(^ + c") 



(24) 



The gaussian path integral over the fields and c" may be easily computed 
with the saddle point method. To this purpose, one has to solve the classical 
equations of motion of these fields: 



c" + LO a d, 



U a C n + iXn 








(25) 

(26) 



5 



In the above equation we have put c° = dc^/dr, = dc^/dr and 

oo a = 2{i\ a + e) (27) 



The solutions of (25-26) satisfying the desired boundary conditions are: 



c n,d( r ) 



a 
An 



1 



a U> a {T-T) 



(28) 
(29) 



Let us note that c" d and c^ cZ do not diverge for large values of r. Moreover, 
it is clear that the boundary conditions are irrelevant in the limit T — > +oo, 
because their contribution vanishes exponentially as e~ 2eT . 
After the field transformation: 



Z r „ becomes: 



Here we have put 



^(r) = < d (r) + < ff (r) 



Z C)n = lim TT f V X %M n , a e iSc n 



(30) 
(31) 



(32) 
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d 



dr 



( a\2 (' v a l 2 



and 



M n , a = J Vc^ ;q V^ q exp i J 



T 



1 -a a , ^ -a -a 



(33) 



(34) 



It is easy to show that the constant factor ftf n>a produced by the integration 
over the "quantum" fields c" ,c^ „ is a just a constant, which is independent 
of u> a and thus can be ignored^. At this point it is possible to perform the 
integration over the pseudo-time r in the action S^ a . The result is: 



qcl _ o;(Xn) T . (Xn) 2 „-2u) i 
°n,a — Li ~ . ,2 V 



' / . 1 * 



Finally, one has to integrate over the variables Xn in %cn' 

7 — 



K + B?) (35) 



VXn exp <^ i 



^7, 



+ (l - e ~W) (5 a + B a, 



(36) 



1 Let us note that this factor is also independent on T. 
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Only the first term in the right hand side of the above equation becomes relevant 
when T becomes very large, as it is in our case. Since 

— < 37 > 

it turns out that, thanks to the presence of the e— term in the action S n of 



Eq. (16), the integrals in T>Xn are convergent. Upon renormalizing the fields 
Xn in Eq. (^) as follows: Xn = Xn^ 1 ^ 2 one finds: 

Zc,n = ]J V^a + it = y / det( J E - H + ie) (38) 

a 

This proves Eq. (]l4| ) as desired. An analogous formula can be derived for ZT 1 . 

Coming back to the original averaging problem, we rewrite Eq. (|9|) in the 
form: 



<G B (q, q'; 0)) - = Jim / V(pP((p)Z c ^Z c ^ 

x J 2^£>#(q)<Kq') exp ^-i J JjL J d D q$(E - ie - H)<j>\ (39) 

The dependence of G^(q, q'; (p) on the disorder fields (p is now explicit and it is 
given by Eq. ( 16 ) , which expresses the inverse of the partition function Z_ as a 
path integral over gaussian fields c n , c n , Xn, fi = 1, 2. At this point it is possible 
to perform the averaging over the disorder fields at least perturbatively. In 
the case of Hamiltonians like that of Eq. (||), the random potential ip\ may be 
integrated out from the partition function with the help of a gaussian integral. 

Let us note that in Eq. ( |39| ) the limit for T — > +oo has been permuted with 
the integration over the disorder fields (p. In a similar way, in the method of 
replicas one needs to permute the limit of vanishing replica's and the average 
with respect to the disorder. The difference between the two approaches is that 
in the present case the limit T — > +oo does not require a complex analytical 
continuation as in the method of replicas and it is mathematically under control. 
In fact, the presence of the variable T is only limited to the partition functions 
Z c ^ n of the bosonic fields c n ,c n ,Xn- These partition functions may always be 
rewritten as in Eq. (|36l), i. e. in terms of standard integrals over the real 
variables x"; which are convergent in the limit T — > +co due to the presence 
of the e term. As a consequence, the permutation of the symbol lini7^ +00 with 
the integrals which are necessary to compute Z c n is allowed. 



4 Conclusions 

Concluding, in this work it has been presented a bosonization approach to dis- 
ordered systems based on the Fermi-Bose principle of which is alternative 
to the supersymmetric method and to the method of replicas. The bosonization 
approach does not contain fermionic degrees of freedom because it involves only 
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bosonic fields. As a consequence, it may be useful in the investigations of sys- 
tems with spontaneous symmetry breaking, which are sometimes complicated 
by the presence of anticommuting variables. 

Let us note that in the case of quantum chromodynamics the Fermi-Bose 
principle used to express the Faddeev-Popov determinant without the help of 
ghost fields leads to a new symmetry, which replaces the usual BRST symmetry. 
It would be thus interesting to check if also the disordered path integral of 
Eq. (|39| ) enjoys an analogous new symmetry, which would replace the fermion- 
boson symmetry of the supersymmetric method. 
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